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Abstract 

Four-quark states may exist as colorless meson-meson molecules or compact systems with two- 
body colored components. We derive an analytical procedure to expand an arbitrary four-quark 
wave function in terms of nonorthogonal color singlet-singlet vectors. Using this expansion we de- 
velop the necessary formalism to evaluate the probability of physical components with an arbitrary 
four-quark wave function. Its application to characterize bound and unbound four-quark states as 
meson-meson, molecular or compact systems is discussed 

PACS numbers: 14.40.Nd,14.40.Lb,14.40.-n 
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I. INTRODUCTION 



The physics of heavy quarks has become one of the best laboratories exposing the hm- 
itations and challenges of the naive quark model and also giving hints into a more mature 
description of hadron spectroscopy. More than thirty years after the so-called November 
revolution heavy meson spectroscopy is being again a challenge. Its formerly comfort- 
able world is being severely tested by new experiments reporting states that do not fit into a 
simple quark- ant iquark configuration It seems nowadays unavoidable to resort to higher 
order Fock space components to tame the bewildering landscape arising with these new 
findings. Four-quark components, either pure or mixed with qq states, constitute a natural 
explanation for the proliferation of new meson states They would also account for the 
possible existence of exotic mesons as could be stable ccnn states, the topic for discussion 
since the early 80 's [4]. 

Four-quark systems present a richer color structure than standard baryons or mesons. 
While the color wave function for standard mesons and baryons leads to a single vector, 
working with four-quark states there are different vectors driving to a singlet color state 
out of colorless or colored quark-antiquark two-body components. Thus, dealing with four- 
quark states an important question is whether we are in front of a colorless meson-meson 
molecule or a compact state, i.e., a system with two-body colored components. While the 
first structure would be natural in the naive quark model, the second one would open a new 
area on the hadron spectroscopy. 

In this manuscript we derive the necessary formalism to evaluate the probability of phys- 
ical channels (singlet-singlet color states) in an arbitrary four-quark wave function. For 
this purpose one needs to expand any hidden-color vector of the four-quark state color 
basis, i.e., vectors with non-singlet internal color couplings, in terms of singlet-singlet color 
vectors. We will see that given a general four-quark state [gi 92^3^4] the above procedure re- 
quires to mix terms from two different couplings, [{qiqs) {q2<l4)] and [(q'i^4)(g2^3)]- If (q'15?2) 
and (ga, ^4) are identical quarks and antiquarks then, a general four-quark wave function 
can be expanded in terms of color singlet-singlet nonorthogonal vectors and therefore the 
determination of the probability of physical channels becomes cumbersome. A particular 
case has been discussed in the literature trying to understand the light scalar mesons as KK 
molecules jsl. This problem has also been found in other fields, as for example in molecular 
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systems [6j. 

The manuscript is organized as follows. In Sec. [TT] the formalism to expand the four- 
quark wave function in terms of singlet-singlet color vectors is derived. Without loss of 
generality, the evaluation of the probabilities is exemplified with the QQnn system. In 
Sec. mil we discuss some examples of bound and unbound states in the charm and bottom 
sectors. Finally, we summarize in Sec. [IV] our conclusions. 



II. FORMALISM 



Given an arbitrary state made up of two quarks, qi and g2, and two antiquarks, 
and g4, its most general wave function will be the direct product of vectors from the color, 
spin, flavor and radial subspaces. We start discussing the color substructure. 



A. Color substructure. 



There are three different ways of coupling two quarks and two antiquarks into a colorless 
state: 

[(gig2)(M4)] = {|3i2334), I612634)} = {|33)i2, 166)^2} (la) 

[(gig3)(g2g4)] = {II13I24), I813824)} = {|ll)c, |88) J (lb) 

[{qiqA){q2q;)] = {II14I23), I814823)} = {|l'l')c, |8'8%} , (Ic) 

being the three of them orthonormal basis. Each coupling scheme allows to define a color 
basis where the four-quark problem can be solved. The first basis, Eq. (fTaj) . being the 
most suitable one to deal with the Pauli principle is made entirely of vectors containing 
hidden-color components. The other two, Eqs. (jlbp and (fTcIl . are hybrid basis containing 
singlet-singlet (physical) and octet-octet (hidden-color) vectors. 

In order to express a four-quark state in terms only of physical components it is necessary 
to define the antiunitary transformation connecting the basis ( ]lb|) and ( fTcl) [3] 



11), = cos a I l'l'), + sina | 8'8'), 

88), = sin a I I'l'), - cos a | 8'8'), , (2) 
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and the projectors on the different vectors: 



P 

Q 



88), , (88 



(3) 



and 



P = |i'i\,(n'| 

Q = I 8'8'), c (8'8' I . 
With these definitions any arbitrary state can be written as 

1^) = P|*) + Q|*) =P\^) + Q\^). 



(4) 



(5) 



One can extract the singlet-singlet components from the octet -octet one, Ql"^), by inserting 
identities l. = P + Q = P + Qm the following iterative manner: 



1^) 



p|^) + pQ|^) + gg|^) = 

P\^) + PQ\^) + PQQ\^) + QQQ\^) = 

P\^) + PQ\^) + PQQ\^) + PQQQl^) + QQQQl^) 



P 



t + QQ + QQQQ 



Q + QQQ + QQQQQ 



(6) 



From the definition of the projectors in Eqs. ([3]) and (jlj) and the antiunitary transformation 
in Eq. one can see that 



QQQ = I 88), , (88 | 8'8'), , (8'8' | 88), , (88 



Therefore, Eq. ([6]) can be rewritten as. 



5'o'\ |2 



)J'g = cos'aQ. (7) 



1^) = P|^) + PQQ [1 + cos^ a + cos^ a + . . .] |^) 
+ PQ[l + cos^ a + cos^ a + ...] |^) . 



Noting that. 



Eq. ([8]) becomes. 



fc=0 



cos a 



2k 



1 



1 — cos^ a 



(8) 



(9) 



I VI/) 



P + PQQ 



1 — COS^ Q! 



1^) 



pg 



1 — cos^ a 



(10) 



This equation can be simplified by considering 



PQP = I 11)^ , (11 I 8'8')^ e (8'8' I 11)^ , (11 I = sin^ a P 



(11) 



and thus 

P + PQQ, 



cos^ a 



P + PQ{l-P) 
PQz 



1 



1 — cos^ a 
P - PQP- 



(12) 



PQ- 



1 



1 — cos^ a 1 — cos^ a 1 — cos^ a 

Therefore Eq. (ITU]) can be finally written in a compact form as 

1 „A,_ 1 



1 — cos^ a 1 — cos^ a 



(13) 



where Pj^^^ and are nonhermitian projection operators on the corresponding singlet- 

singlet subspaces (see Appendix |X] for proof of their properties). This expression demon- 
strates that any octet-octet color component can be expanded, in general, as an infinite 
sum of singlet-singlet color states Q]. 

To obtain hermitian operators one can repeat the same procedure using the projectors 
on I 1'1% and | 8'8% given in Eq. (HD, 



1^) 



Q\^) + QP\^) + PP\^) = 

Q\^) + QP\^) + QPP\^) + PPP\^) = 

Q\^) + QP\^) + QPP\^) + QPPP\^) + PPP\^) 



Q 



1 + PP + pppp + ... 



p + ppp + ppppp 



(14) 



where 



PPP = I I'l')^ , (I'l' I 11)^ , (11 I I'l'), c (I'l 

what allows to rewrite Eq. (fT4|) 



(I'l' I 11)J^ P = cos^a P, (15) 



1^) = Ql^) + gpp [l + cos2 a + cos^ a 
+ QP [l + cos^ a + cos^ « + ...] 1^) 



...]\^) 



(16) 



Using Eq. ([9]) one can finally write 



I*) 



Q + QPP 



1 — cos^ a 



QP 



1 — cos^ a 



(17) 
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This equation can be simplified noting that 



QPQ = I 8'8'), c (8'8' I 11)^ c (11 I 8'8% , (8'8' | = sin^ a Q , 



then 



Q + QPP- ^ = Q + QP{1-Q)- ^ 

1 — cos^ a 1 — cos^ a 



(18) 



(19) 



= QP 

arriving to the compact notation 

1 



1 — cos^ a 



Q-QPQ 



1 — cos^ a 



QP 



1 — cos^ a 



1 — cos^ a 1 — cos^ a 



,^QP\^) = v^,^^l\^) + v^,!',?^m 



(20) 



where Pj^^^^ and P^^!(^ are nonhermitian projection operators on the corresponding singlet- 
singlet subspaces. 

Combining Eqs. f|T3|l and fl20l) one can write any arbitrary state in the following form, 



1^) 



+ 



PQ 
QP 



1 — cos^ a 
1 

1 — cos^ a 



PQ 
QP 



1 



1 — cos^ a 
1 

1 — cos^ a 



(21) 



or equivalently 



|v^) = UpQ + Qp) L^lvl/) 

2 \ /I — cos"^ a 

+ Upq + qp)- — ^1^). 

2 \ /I — cos"^ a 



(22) 



Thus, one arrives to two hermitian operators that are well-defined projectors on the two 
physical singlet-singlet color states 

1 



PQ + QP 



PQ + QP 



2(1 — cos^ a) 
1 



(23) 



2(1 — cos^ a) 
and finally, 

l^)=^|ii)J^)+^|i'i')J*)- (24) 

Thus, given an arbitrary state its projection on a particular subspace E is given by 
= -PeI^)- Thus, the probability of finding such an state on this subspace is 



(25) 
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(26) 



Therefore, once the projection operators have been constructed, Eq. (l23l) . the probabihties 
for finding singlet-singlet components are given by, 

Pl^>([ll]) = (^|P|n)J^> 
pl*)([l'l']) = I V\vv)^ I . 

Using Eq. (!23|) one arrives to 

^''^M)- 2(l-cos^.) + 
Pl^)([l'l']) = ^(^3^^ [(^ I PQ I ^) + I QP I ^)] (27) 

where it can be easily checked that pl*>([ll]) + pl'^>([l'l']) = 1. 



B. Spin substructure 

For a four-quark state one has three different total spins: 0, 1 and 2. The St = 2 case is 
trivial, because the basis is one-dimensional. Let us discuss the other two possibilities. For 
St = the spin basis, in analogy with Eqs. (1), are given by: 

[(^i^2)5..(^3^4)53jo ^{miMim 

and the corresponding spin projectors 

Ps = |00),,(00| 

Qs = |ii).s(ii| 

P, = |0'0%,(0'0'| 

It is important to note that the projectors used in the color space determine the coupling 
in the spin space. Thus, introducing the corresponding spin projectors in Eq. 
arrives to 



(28a) 
(28b) 
(28c) 

(29) 



one 



\^) = P|n>^ (Ps + Qs) m + P|i'i'>. [Ps + Qs) m = 

= Viu)^Ps\^) +Vin)^Qs\'^) +Virv)/s\'^) +Virv)^Qs\^) = 

= 'P|11)^,|00>J*) +P|ll)^,|ll>J*) +^|1'1')„|0'0')J*) +'^|i'i')c-|i'i'>J*) 



(30) 



where ^111)^,100)8 s^nd 'P|i'i')^,|o'o')s stand for the projectors on the physical state made up of 
two S = qq mesons, and "^111)^,111)^ and 'P|i'i')^,|i'i')s ^^e projectors on the physical state 
made up of two S = 1 qq mesons. 

Following our discussion in Subsection III Al the probabilities are given by 

P[|11)J00),] = {^\Vin)joo).m (31) 

P[|11)J11).] = {^\Vin)jn).m 

P[|1T)J0'0').] = {^\Viry^jo'o')M 

^[|1'1')J1'1%] = (*l^|i'i%|i'i'>J*), 

and therefore, the total probabilities of finding a physical state made up of two S = qq 
mesons will be given by 

Pmm = P[\ 11), |00),] + P[\ I'l'), |0'0').] (32) 
and correspondingly the total probability of a physical state made up of two S = 1 qq states 

Pm*m*= P[|ll)Jll).]+P[|lT)jn%]. (33) 
In the St = 1 case the spin basis are 

[{SlS,)s,As3S,)s.Jl = {|01)i^ 110)12, 111)12} (34^) 

[{siS3)sAs2S,)sJi = {|01)., |10)„ 111),} (34b) 
[isis,)s,As2S3)s,s]i = {|0'1%, |1'0%, |1'1%} (34c) 

and the corresponding projectors. 



Ps = 


|01).s(01| 


Qs = 


|10)..(10| 


Ws = 


|ll)ss(ll| 


Ps = 


|0'1%.(0'1'| 


Qs = 


|1'0%.(1'0'| 


= 


|i'i%s(i'i'| 



(35) 
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Following the same procedure as in the St = case one arrives to 

1^) = Viu)^ {Ps + Qs + Ws) 1^) + P|ia'>. (A + Qs + Ws) 1^) = (36) 
= P|n) P.lvl^) + + Vin)Ws\^) 

= "PudJodM +Viu)jio)M +Viu)ju)M 

+ ^|1'1'>„|0'1'>J^) + ^|1'1'>„|1'0'>J^) + '^|1'1'>„11'1'>J^) > 

where P|ii)^,|oi),, '^|i'i'>„|o'i'>., ^|ii>„|io>,, and P|i'i'>„|i'o'>. stand for the projectors on the 
physical state made up of one S = and one S = 1 qq mesons and ^111)^,111)^ and Pin)^,!!'!')^ 
for the projectors on the physical state made up of two S = 1 qq mesons. 
Finally, the probabilities can be expressed as 

P[|11)J01).] = mVinuoDJ^) (37) 

P[|11)J10).] = {^\Vinuio).m 

P[|11)J11),] = (vl/|P|n)jn)J*) 

P[|l'l')c|0'l%] = mVirr^jo'i%m 

P[|1'1')J1'0%] = {^\Viyr)ji'o')M 

P[| 1'1')J1'1%] = (M/|P|ia')Ji'i')J^), 

and therefore, the total probability of a physical state made up of one S = and one S = 1 
qq meson will be given by 

Pmm*=P[\ 11)J01).]+P[| 11)J10).]+P[| 1'1')J0T).]+P[| 1Y)J1'0%] (38) 

and the total probability of a physical state made up of two S = 1 qq mesons by 

Pm.m* = P[|11)J11).]+P[|1'1')J1'1%]. (39) 

C. Flavor substructure 

The previous discussion about the color and spin substructure is general and valid for 
any four-quark state. For the flavor part one find several different cases depending on the 
number of light quarks. Although the present formalism can be applied to any four-quark 
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state, it becomes much simpler whether distinguishable quarks are present. This would be, 
for example, the case of the nQnQ system, where the Pauli principle does not apply. In this 
system the basis (llbp and (ITcl) are distinguishable due to the flavor part, they correspond 

to [(iiu)/i4 (^2^3)723]/ = [(?c)i/2(cg)i/2]/ and [{iih) iri{i2H) I2a]i = [(9^)/(cc)o]/, and therefore 
they are orthogonal. This makes that the probabilities can be evaluated in the usual way 
for orthogonal basis as has been done in Ref. 9|. 

Non-orthogonal basis are necessary for the following cases: QQnQ', QQ'nn, Qnnn and 
nnnn {Q may be equal to Q') or their corresponding antiparticles. The isospin basis are: 



QQnQ' 



QQ'nn 



[(^1^2)7,2(^4)734] I =|0|)f (40a) 

[(^1^3)7i3(«2^4)724]i = l|0)/ (40b) 

[(^1^4)7,4 (^2^3)723]! =|0'i')/- (40c) 



[(^1^2)712 (^4)734)7 = |0/)f (41a) 
[(«i«3)7i3 (^2^4)724)7 = III)/ (41b) 
[(^1^4)714 (^2^3)723)7 =l|'|')/- (41c) 



Qnnn 

~ I =1/2 



I = 3/2 



[(^1^2)7,2(^4)734]! ={||0)}M|l)f} (42a) 
[(^3)7.3(^2^4)724]! ={||0)/,||1)/} (42b) 
[(^4)7,4 (^2^3)723]! ={l|'0')/,||'l')/}- (42c) 



[(^1^2)712 (^3^4)734] I = l|l)f (43a) 
[(«i«3)7i3 (^2^4)724] I = l|l)/ (43b) 
[(«i«4)7i4 (^2^3)723] I =l|'l')/- (43c) 
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nnnn 



1 = 

[{^l^2)IA^3^4)I,,]o = {|00)f,|ll)f } (44a) 

[(h^3)/.3(^2^4)/Jo ={|00);,|11)^} (44b) 

[{^l^,)IA^2^3)I,■s]o = {|0'0')/, • (44c) 



[(^1^2)/,.(Z3^4)/3J1 ={|01)}MlO)f,|ll)f} 
[(V3)/i3(^2^4)7,Jl ={|01);,|10);,|11);} 
[(M4)/,,(^2^3)7.3]l ={|0'l')/,|l'0')y,|l'l')^} 



(45a) 
(45b) 
(45c) 



-1 = 2 

[i^l^2)IA^3^A)IJ2 = (46a) 

[(H^3)/i3(«2^4)/24]2 = 111)/ (46b) 
[(n^4)/,,(^2^3)/,3]2 =|1'1')/- (46c) 

For those cases where the basis is one-dimensional the recouphng among the three different 
basis introduced in Eq. ([1]) is straightforward. For those cases where the basis in not one- 
dimensional one should follow the procedure described in Sect. IIIBI 

D. Radial substructure 

In order to derive the probability of the physical channels one has finally to analyze the 
symmetry of the radial wave function. Such analysis will depend on the particular state 
chosen. Without loss of generality we will exemplify the procedure with the particular case 
of the QQnn system. Any other four-quark system discussed in Sect. Ill CI could be analyzed 
in the same manner. Let us start with the St = state, whose most general wave function 
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reads 

\^) = \R,mi'mi'\oi)f + \R,mi'\ii)i'\oi)f (47) 

+ |i?3)|66)J2|00)f|0/)f + |i?4)|66)^2|ll)f|0/)f , 

where |-Ri), 1-^2); l-^s); and I-R4) are radial wave functions that due to symmetry properties 
satisfy {Ri\Ri) + (i?2|i?2) + {RslRs) + {Ri\Ri) = 1, (i?i|i?2) = {RilRs) = {R2\Ri) = 
{RsIRa) = 0, (RiIRa:) ^ 0, and (i?2|^?3) 7^ 0. Applying Eqs. ([32]) and ([33]) one obtains 

Pmm = P[|ll)e|00),]+P[|l'l')c|0'0'),] (48) 
^ '1 + 2(^21^2) + 2(i?4|i?4)) + ^ ((P1IP4) + (P2IP3)) 
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Pa/.a/* = P[|ll)c|ll).]+P[|l'l')c|l'l') 



= i (1 + 2(Pi|Pi) + 2(P3|P3)) - ^ ((P1IP4) + (P2IP3)) . 
4 o 

Finally, the QQnn St = ^ most general wave function reads 

m = \R^ml'\01)l'\0I)f + |P2)|33)f |10)f |0/)f + |P3)|33)^2|ll)^2|0/)f (49) 

where \Ri) are radial wave functions that due to symmetry properties satisfy ^^=i(Pi|Pi) = 
1 and all the cross products are zero except for (P1IP5) and (^21^4)- Applying Eqs. (138|1 
and ( l39l) one gets 

^ ' l + (P3|P3) + (P6|P6)-^((Pl|P5) + (P2|P4)) 1 (50) 



MM" 



2 \ ^ ^ 2 



Pa/.m. = ^ ( ^ " ^^^l^a) - (PeiPe) + ^ ((P1IP5) + (P2IP4)) 



III. SOME ILLUSTRATIVE EXAMPLES 



In the previous section we have derived the analytical expansion of an arbitrary four- 
quark state wave function in terms of a non-orthogonal basis containing only physical chan- 
nels. The calculation of the probabilities has been exemplified with the QQnn system. 

We now apply this formalism to discuss the four-quark nature: unbound, molecular or 
compact states, of some illustrative examples. The same discussion could be done with any 
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FIG. 1: -ff-type Jacobi vectors. 1,2 stand for heavy quarks and 3,4 for light antiquarks. 




other four-quark state just changing the couphng in isospin space, see Subsect. Ill CI The 
results we are going to discuss have been obtained solving the four-body problem by means 
of a variational method using as trial wave function the most general linear combination of 
gaussians [l^ . The accuracy of the variational approach has been tested by comparing with 



results obtained by means of the hyperspherical harmonic expansion 



ll| . Both approaches 



are i n g ood agreement The interaction between the quarks is taken from the model of 



Ref . [l^ . The same interacting potential used to calculate the four-quark energy is used to 
calculate the mass of the thresholds, i.e., the meson masses. 

The stability of a four-quark state can be analyzed in terms of A^;, the energy difference 
between its mass and that of the lowest two-meson threshold, 

A£; = E4,-E(Mi,M2), (51) 

where E^g stands for the four-quark energy and E{Mi, M2) for the energy of the two-meson 
threshold. Thus, A^; < indicates all fall-apart decays are forbidden, and therefore one has 
a proper bound state. A^; > will indicate that the four-quark solution corresponds to an 
unbound threshold (two free mesons). Thus, an energy above the threshold would simply 
mean that the system is unbound within our variational approximation, suggesting that the 
minimum of the Hamiltonian is at the two-meson threshold. Another helpful tool analyzing 
the structure of a four-quark state is the value of the root mean square radii: (x^)^/^, (y^)^/^, 
and (2;^)^/^. They correspond to the Jacobi coordinates given in Fig{Tl Compact four-quark 
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TABLE I: Probabilities for ccnn states with quantum numbers {St, I) = (1,1) and (1,0). The 
notation used stands for P[|{Color state}) |{Spin state})] where {Color state} corresponds to the 
basis vectors given in Eqs. ([T]) and {Spin state} to the ones given in Eqs. (|34p . Flavor component 



\OI)f, 1^^)/, and 1^'^')/ is understood. We list in Appendix [Cl a summary of the expressions used. 



(1,1) (1,0) 


(1,1) (1,0) 


(1,1) (1,0) 


p 


|33)i2|01),] 


0.000 0.875 


P 


1 11 


)e 


01)s] 


0.277 0.094 


P[ 


1 


O'l') 


s] 


0.277 0.094 


p 


|33)i2|10).] 


0.000 0.006 


P 




)e 


10)s] 


0.277 0.094 


P[ 


1 


I'O'] 


s] 


0.277 0.094 


p 


|33)^2|11),] 


0.333 0.000 


P 


1 11 


>c 


ll)s] 


0.002 0.186 


P[ 


1 


I'l'] 


s] 


0.002 0.186 


p 


|66)i2|01),] 


0.000 0.090 


P 


1 88 


>c 


01)s] 


0.222 0.156 


P[ 


1 8'8'), 


O'l'] 


s] 


0.222 0.156 


p 




0.000 0.029 


P 


1 88 


>e 


10)s] 


0.222 0.156 


P[ 


1 8'8'), 


I'O'] 


s] 


0.222 0.156 


p 


|66)^2|n),] 


0.667 0.000 


P 


1 88 


)e 


ll)s] 


0.000 0.314 


P[ 


1 8'8'), 


I'l'] 


s] 


0.000 0.314 




pm?] 


0.333 0.881 




^[| 


11 


)J 


0.556 0.374 




P[|l'l 


')c] 




0.556 0.374 




p[\m)f\ 


0.667 0.119 




^[| 


88 


)J 


0.444 0.626 




P[ \ 8'8 


')J 




0.444 0.626 



states can be distinguished from two free mesons by means of their root mean square radius 



RMS 



4<?(2q) 



1/2 



(52) 



and in particular, their corresponding ratio, 



A 



RMS, 



R 



4q 



(53) 



RMSa/i + RMSm2 

where RMSmi + RMSa./2 stands for the sum of the radii of the mesons corresponding to the 
lowest threshold. 



We show in Table [T] the probabilities in color and spin space obtained for two ccnn states. 
The first one, with quantum numbers {St, I) = (1,1), is unbound, while the second one, 
{St, I) = (1,0), is bound. We give the probabilities in the three different rearrangements 
in color space defined in Eqs. ([T]). Let us note that in the three color rearrangements of 
Eqs. ([T]), the hidden-color vectors (|33), |66), |88), and |8'8')) contain probability of physical 
channels as we have discussed in Sect. Ull It is possible to prove from simple group theory 
arguments that for a system composed of two identical quarks {QQ) and two identical 
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TABLE II: Four-quark state properties for selected quantum numbers. All states have positive 
parity and total orbital angular momentum L = 0. Energies are given in MeV and distances in fm. 
The notation M1M2 \i stands for mesons Mi and M2 with a relative orbital angular momentum 
I. P[|33)J^(|66)J^)] stands for the probability of the 33(66) components given in Eq. (fTa|) and 
P[\ 88)^)] for the 11(88) components given in Eq. (|lb[) . Pmm, Pmm*, and Pm*m* have 

been calculated as explained in text. 



{St, I) 


(0,1) 


(1,1) 


(1,0) 


(1,0) 


(0,0) 


Flavor 


ccnn 


ccnn 


ccnn 


bbnn 


bbnn 


Energy 


3877 


3952 


3861 


10395 


10948 


Threshold 


DD \s 


DD* \s 


DD* \s 


BB* \s 


BiB \p 


Aij 


+5 


+15 


-76 


-217 


-153 


P[|33)i2] 


0.333 


0.333 


0.881 


0.974 


0.981 


P[|66)i2] 


0.667 


0.667 


0.119 


0.026 


0.019 


LI / CJ 


0.556 


0.556 


0.374 


0.342 


0.340 


^'[l 88)J 


0.444 


0.444 


0.626 


0.658 


0.660 


Pmm 


1.000 








0.254 


Pmm* 




1.000 


0.505 


0.531 




Pm*m* 


0.000 


0.000 


0.495 


0.469 


0.746 


(x2)l/2 


60.988 


13.804 


0.787 


0.684 


0.740 




60.988 


13.687 


0.590 


0.336 


0.542 


(z2)l/2 


0.433 


0.617 


0.515 


0.503 


0.763 


RM S^q 


30.492 


6.856 


0.363 


0.217 


0.330 




69.300 


11.640 


0.799 


0.700 


0.885 



antiquarks {nn), the octet-octet component probability of the wave function either in the 
(]lb|) or (JTcj) arrangements is restricted to the interval [1/3,2/3], see Appendix [Bl Do these 
numerical and analytical results prove the unavoidable presence of important hidden-color 
components in all QQnn states regardless of their binding energy? We shall see that the 
answer is no. 
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To respond this question we summarize in Table [TTl the results obtained for several different 
four-quark states, among them those used in Table [H making use of the formal development 
of Sect, [m As shown in Table IIT| independently of their binding energy, all of them present a 
sizable octet-octet component when the wave function is expressed in the fllb|) coupling. Let 
us first of all concentrate on the two unbound states, A^; > 0, one with St = and one with 
St = 1, given in Table HTl The octet-octet component of basis (lb) can be expanded in terms 
of the vectors of basis (Ic) as explained in the previous section. Thus, once expressions (42) 
and (44) are considered one finds that the probabilities are concentrated into a single physical 
channel, Pmm or Pmm*- In other words, the octet-octet component of the basis (lb) or (Ic) 
is a consequence of having identical quarks and antiquarks. Thus, four-quark unbound states 
are represented by two isolated mesons. This conclusion is strengthened when studying the 
root mean square radii, leading to a picture where the two quarks and the two antiquarks 
are far away, (x^)^/^ ^ 1 fm and (y^)^^^ ^ 1 fm, while the quark-antiquark pairs are located 
at a typical distance for a meson, (z^)^/^ < 1 fm. 

Let us now turn to the bound states shown in Table HTl A^; < 0, one in the charm sector 
and two in the bottom one. Contrary to the results obtained for unbound states, when 
the octet-octet component of basis (lb) is expanded in terms of the vectors of basis (Ic), 
equations (42) and (44) drive to a picture where the probabilities in all allowed physical 
channels are relevant. It is clear that the bound state must be generated by an interaction 
that it is not present in the asymptotic channel, sequestering probability from a single 
singlet-singlet color vector due to the interaction between color octets. Such systems are 
clear examples of compact four-quark states, in other words, they cannot be expressed in 
terms of a single physical channel. Moreover, as can be seen in Table [TTl their typical sizes 
point to compact objects 20% smaller than a standard two-meson system. 

We have studied the dependence of the probability of a physical channel on the binding 
energy. For this purpose we have considered the simplest system from the numerical point of 
view, the {St, I) = (0, 1) ccnn state. Unfortunately, this state is unbound for any reasonable 
set of parameters. Therefore, we bind it by multiplying the interaction between the light 
quarks by a fudge factor. Such a modification does not affect the two-meson threshold while 
it decreases the mass of the four-quark state. The results are illustrated in Fig. [2] (Pmm) 
and Fig.[3](A/j, {x^^^, and {z^y/"^). In Fig.[2]it is shown how in the A^ ^ limit, 

the four-quark wave function is almost a pure single physical channel. Close to this limit 
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FIG. 2: Pmm as a function of A^;. 
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one would find what could be defined as molecular states. In Fig. [3] we see how the size of 
the four-quark state increases when A^; — > 0. It can be observed that when the probability 
concentrates into a single physical channel {Pmm ~^ 1) the size of the system gets larger 
than the sum of two isolated mesons (Fig. [3] left panel). In Fig. [3] (right panel) we identify 
the subsystems responsible for increasing the size of the four-quark state. Quark-quark 
((x^)^/^) and antiquark-antiquark ((y^)^^^) distances grow rapidly while the quark-ant iquark 
distance ((z^)^/^) remains almost constant. This reinforces our previous result, pointing to 
the appearance of two meson like structures whenever the binding energy goes to zero. This 
illustrative example emphasizes the importance of performing a simultaneous analysis both 
of energy and wave function in order to detect bound states in the vicinity of a two-meson 
threshold. 

To illustrate in more detail the differences observed in the calculated four-quark wave 
functions we depict in Fig. |4]the position distributions defined as 



where Ri{ra^rp,r^) are the radial wave functions introduced in Eqs. ( 1471) and ( 149|) . We 
present results for an unbound, a molecular and a bound state, showing the position distri- 
bution for the different planes (rQ,,r^) = (2, x), (2;, ?/), and {x,y). Clear differences among 
them can be observed. The position distribution for the unbound case spreads in the x and 




(54) 
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FIG. 3: Ar (a) and (x2)V2(solid line), (y2)i/2 (dashed line), and (z2)i/2(dashed-dotted line) (b) 
as a function of A^. 
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y coordinates up to 60 fm, while the bound and molecular systems are restricted to the 
region below 3 fm (molecular) and 1 fm (bound). In the [x^y) plane the unbound state is 
so widely spread that the values for the position distribution are three orders of magnitude 
lower than in the (-z, y) and (2,x) cases, and therefore they will not appear in the picture 
unless artificially magnified. In the case of the molecular state a long range tail propagating 
in the x = y region can be observed contrary to the constrained values obtained for bound 
systems. 

The conclusions derived are independent of the quark-quark interaction used. They 
mainly rely on using the same hamiltonian to describe tensors of different order, two and four- 
quark components in the present case. Dealing with a complete basis, any four-quark bound 
deeply bound state has to be compact. Only slightly bound systems could be considered as 
molecular. Unbound states correspond to a two-meson system. A similar situation would be 
found in the two baryon system, the deuteron could be considered as a molecular like state 
with a small percentage of its wave function on the AA channel, while the if— dibaryon 
would be a compact six-quark state. Working with central forces, the only way of getting a 
bound system is to have a strong interaction between the constituents that are far apart in 
the asymptotic limit, quarks or antiquarks in the present case. In this case the short-range 
interaction will capture part of the probability of a two-meson threshold to form a bound 
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FIG. 4: Position distribution corresponding to unbound {St = 1, / = 1, ccnn), molecular (St = 0, 
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state. This can be reinterpreted as an infinite sum over physical states. This is why the 
analysis performed in the present manuscript is so important before concluding the existence 
of compact four-quark states beyond simple molecular structures. 

If the prescription of using the same hamiltonian to describe all tensors in the Fock space 
is relaxed, new scenarios may app ear. Among them, the inclusion of many-body forces is 



particularly relevant. In Ref. 



13| the stability of QQnn and QQnn systems was analyzed 



in a simple string model considering only a multiquark confining interaction given by the 
minimum of a fiip-fiop or a butterfly potential in an attempt to discern whether confining 



19 



interactions not factorizable as two-body potentials would influence the stability of four- 
quark states. The ground state of systems made of two quarks and two antiquarks of equal 
masses was found to be below the dissociation threshold. While for the cryptoexotic QQnn 
the binding decreases when increasing the mass ratio mQ/nin, for the flavor exotic QQnn 



the effect of mass symmetry breaking is opposite. Ot 
many-body forces, like many-body color interactions 
body interactions [l^, are considered. 



lers scenarios may emerge if different 



141 ] or 't Hooft instanton-based three- 



IV. CONCLUSIONS 



In this work we have developed the necessary formalism to express the wave function of 
a general four-quark state in terms of physical channels, i.e., those constructed by using 
color singlet states. Once this is done the four-quark wave function is expressed in terms 
of nonorthogonal vectors and hence the traditional way to compute probabilities needs to 
be generalized. We have obtained expressions to evaluate such probabilities for all possible 
nontrivial four-quark states containing two heavy antiquarks and two light quarks. We have 
applied these expressions to illustrative cases, where the difference among unbound, compact 
and molecular four-quark states has been made evident. The importance of performing a 
complete analysis of the system, energy and wave function, in the vicinity of a two-meson 
threshold has been emphasized. 
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APPENDIX A: PROJECTORS PROPERTIES 



The following properties are proved for one particular set of projectors, Pj^^^ and V^J^}-^ . 
The same procedure can be followed in all the remaining cases. By construction, see Eq. 
(24), they span the complete space, V^.^x^ + = 1- We demonstrate that we have 
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(Al) 

88)^e(88|l'l%,(l'l'|88)^e(88| 
88)^ e (88 I U(1'1'|88)J2 
88)^ e(88 I sin^a 



APPENDIX B: MINIMUM AND MAXIMUM VALUE FOR THE OCTET 
OCTET COMPONENT PROBABILITY. 

Without loss of generality we consider the St = case. The Pauli Principle requires that 
the radial wave functions \Ri) in Eq. (41) have well-defined permutation properties under 
the exchange of quarks and that of antiquarks, i.e., symmetric (S) or antisymmetric (A). 
Hence, \Ri) must be antisymmetric under the exchange of the identical quarks and also 
under the exchange of antiquarks what we will denote by \Ri{AA)). Similarly for the other 
components: \R2{SS)), iR-s^SS)), and \R4^{AA)). The transformations from (fTal) to ( flbll 
and from (l28aj) to (l28bl) are 

= 7f iii)^-vf '^^^^ (^^^ 

and 

= ^|00). + ^|11). (B2) 
= ^100). -^|11). 



constructed idempotent operators, 

l'l')cc(l'l' 



1 



1 — cos^ a 



\l — cos^ a J 



1 — cos^ a 



'|l'l')cc(l'l' 



1 — cos"^ a I 
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and therefore Eq. (41) can be written as 

\^) = |ll)jOO),|^i);^{|i?i) + v^|i?2) + v^|i?3) + v^|i?4)}+ (B3) 

I ii)Jii).|^^)/^ {v^|i?i) - |i?2) + Vem - V2\R,)} + 

I 88), I00)s|^^)/^ {-V2|i?i) - V6\R2) + \R3) + V3\R,)] + 

I 88)e {-V6\Ri) + v^|i?2) + VslRs) - \R,)} . 

Considering the symmetry properties of the radial part of the wave function, the probabihties 
are calculated as 

P[\ 11), |00),] = ^ {{Ri\Ri) + 3(i?2|i?2) + 2{Rs\Rs) + QiR^lR^) (B4) 

+ 2^6 Re ((i?i|i?4) + {R2\R3))} 
P[|11)J11).] = ^{3{R,\R,) + {R2\R2) +Q{R3\R3) + '2{R4R4) 

- 2V6Re ({i?i|i?4) + (i?2|i?3))} 

P[|88)J00),] = ^{2(i?i|i?i)+6(i?2|i?2) + (i?3|i?3)+3(i?4|i?4) 

- 2^6 Re ((i?i|i?4) + (i?2|i?3))} 
P[|88)J11),] = ^{6{Ri\Ri) +2{R2\R2) + ^R3\R3) + {RaW 

+ 2^6 Re ((P1IP4) + (P2IP3))} • 

Thus, 

^[|11)J = P[|11)J00).]+P[|11)J11),] (B5) 

= ^ {4(Pi|Pi) + 4(P2|P2) + 8(P3|P3) + 8(P4|P4)} 

P[|88)J = P[|88)J00).]+P[|88)J11).] 

= ^ {8(Pi|Pi) + 8(P2|P2) + 4(P3|^3) + 4(P4|P4)} . 

By construction P[|33),^2] ^ {R^\R^) + (P2IP2) and P[|66)J2] = (P3IP3) + (P4IP4) with 
P[|33)J^] + P[|66)J^] = 1. Therefore Eqs. fIBSP can be expressed as 

^[|88)J = 1{2-P[ml']} . 

and since P[|66)J^] G [0, 1] is normalized, a minimum (1/3) and a maximum (2/3) value for 
P[| 11) J and P[| 88) J do exist. 
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APPENDIX C: PROBABILITIES FOR DIFFERENT CHOICES OF BASIS. 



P[ 
P[ 

P[ 
P\ 
P\ 
P[ 
P[ 
P[ 



The 5't = case is given in Eqs. (IBSp . for the sake of completeness note that 
11)J00)J = P[| 1T)J0'0%], P[|11)J11).] = P[\VV)^l'l')sl P[|88)J00).] = 
8'8')J0'0%],P[| 88)J11),] = P[\ rr)Jl'l%]. For St = 1 one has 

11)J01).] = P[|1'1')J0'1'),]= (CI) 
^ (^1 - \{Ri\Ri) - \{R2\R2) + {R^\R^) - V2[(Pi|P5) + {R2\Ra)] 

11)J10).] = P[|1'1')J1'0%] = 

^ (^1 - ^(PilPi) - \{R2\R2) + {R,\R,) - \/2[(Pl|P5) + {R2\Ra)]^ 

11)J11).] = ^[|l'l')e|l'l').] = 

^ ((PllPl) + (i?2|i?2) + 2(i?4|i?4) + 2(i?5|i?5) + 2v^[(Pi|P5) + (P2IP4)] 



^)J01).] = P[|8'8')J0T).] 
^)J10),] = P[|8'8')J1'0%] 



1 - ]^{RARa) - \{R,\R,) + (i?3|/23) + V2\{R,\R,) + (P2|/24)] 



^ - ^(P4|P4) - \{R^\R^) + (i?3|i?3) + V2[{R^\R,) + (P2|i?4)] 



^)J11),] = P[|8'8')J1Y),] 



^ (2(Pi|Pi) + 2(P2|i?2) + {Ra\Ra) + (i?5|i?5) - 2y2[(Pi|P5) + (P2IP4)] 
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